Disentangling top-down drivers of mortality underlying diel population dynamics of Prochlorococcus in the North Pacific Subtropical Gyre.

Photosynthesis fuels primary production at the base of marine food webs. Yet, in many surface ocean ecosystems, diel-driven primary production is tightly coupled to daily loss. This tight coupling raises the question: which top-down drivers predominate in maintaining persistently stable picocyanobacterial populations over longer time scales? Motivated by high-frequency surface water measurements taken in the North Pacific Subtropical Gyre (NPSG), we developed multitrophic models to investigate bottom-up and top-down mechanisms underlying the balanced control of Prochlorococcus populations. We find that incorporating photosynthetic growth with viral- and predator-induced mortality is sufficient to recapitulate daily oscillations of Prochlorococcus abundances with baseline community abundances. In doing so, we infer that grazers in this environment function as the predominant top-down factor despite high standing viral particle densities. The model-data fits also reveal the ecological relevance of light-dependent viral traits and non-canonical factors to cellular loss. Finally, we leverage sensitivity analyses to demonstrate how variation in life history traits across distinct oceanic contexts, including variation in viral adsorption and grazer clearance rates, can transform the quantitative and even qualitative importance of top-down controls in shaping Prochlorococcus population dynamics.

Widespread yet persistent low abundance of TIM5-like cyanophages in the oceans.

Ocean ecosystems are inhabited by a diverse set of viruses that impact microbial mortality and evolution. However, the distribution and abundances of specific viral lineages, particularly those from the large bank of rare viruses, remains largely unknown. Here, we assessed the diversity and abundance of the TIM5-like cyanophages. The sequencing of three new TIM5-like cyanophage genomes and environmental amplicons of a signature gene from the Red Sea revealed highly conserved gene content and sequence similarity. We adapted the polony method, a solid-phase polymerase chain reaction assay, to quantify TIM5-like cyanophages during three 2000 km expeditions in the Pacific Ocean and four annual cycles in the Red Sea. TIM5-like cyanophages were widespread, detected at all latitudes and seasons surveyed throughout the photic zone. Yet they were generally rare, ranging between <100 and 4000 viruses·ml-1 . Occasional peaks in abundance of 10- to 100-fold were observed, reaching 71,000 viruses·ml-1 . These peaks were ephemeral and seasonally variable in the Red Sea. Infection levels, quantified during one such peak, were very low. These characteristics of low diversity and abundance, as well as variable outbreaks, distinguishes the TIM5-like lineage from other major cyanophage lineages and illuminates that rare virus lineages can be persistent and widespread in the oceans.

Cyanophages from a less virulent clade dominate over their sister clade in global oceans.

Environmental virus communities are highly diverse. However, the infection physiology underlying the evolution of diverse phage lineages and their ecological consequences are largely unknown. T7-like cyanophages are abundant in nature and infect the marine unicellular cyanobacteria, Synechococcus and Prochlorococcus, important primary producers in the oceans. Viruses belonging to this genus are divided into two distinct phylogenetic clades: clade A and clade B. These viruses have narrow host-ranges with clade A phages primarily infecting Synechococcus genotypes, while clade B phages are more diverse and can infect either Synechococcus or Prochlorococcus genotypes. Here we investigated infection properties (life history traits) and environmental abundances of these two clades of T7-like cyanophages. We show that clade A cyanophages have more rapid infection dynamics, larger burst sizes and greater virulence than clade B cyanophages. However, clade B cyanophages were at least 10-fold more abundant in all seasons, and infected more cyanobacteria, than clade A cyanophages in the Red Sea. Models predicted that steady-state cyanophage abundances, infection frequency, and virus-induced mortality, peak at intermediate virulence values. Our findings indicate that differences in infection properties are reflected in virus phylogeny at the clade level. They further indicate that infection properties, together with differences in subclade diversity and host repertoire, have important ecological consequences with the less aggressive, more diverse virus clade having greater ecological impacts.

Viruses affect picocyanobacterial abundance and biogeography in the North Pacific Ocean.

The photosynthetic picocyanobacteria Prochlorococcus and Synechococcus are models for dissecting how ecological niches are defined by environmental conditions, but how interactions with bacteriophages affect picocyanobacterial biogeography in open ocean biomes has rarely been assessed. We applied single-virus and single-cell infection approaches to quantify cyanophage abundance and infected picocyanobacteria in 87 surface water samples from five transects that traversed approximately 2,200 km in the North Pacific Ocean on three cruises, with a duration of 2-4 weeks, between 2015 and 2017. We detected a 550-km-wide hotspot of cyanophages and virus-infected picocyanobacteria in the transition zone between the North Pacific Subtropical and Subpolar gyres that was present in each transect. Notably, the hotspot occurred at a consistent temperature and displayed distinct cyanophage-lineage composition on all transects. On two of these transects, the levels of infection in the hotspot were estimated to be sufficient to substantially limit the geographical range of Prochlorococcus. Coincident with the detection of high levels of virally infected picocyanobacteria, we measured an increase of 10-100-fold in the Synechococcus populations in samples that are usually dominated by Prochlorococcus. We developed a multiple regression model of cyanophages, temperature and chlorophyll concentrations that inferred that the hotspot extended across the North Pacific Ocean, creating a biological boundary between gyres, with the potential to release organic matter comparable to that of the sevenfold-larger North Pacific Subtropical Gyre. Our results highlight the probable impact of viruses on large-scale phytoplankton biogeography and biogeochemistry in distinct regions of the oceans.

Genetic engineering of marine cyanophages reveals integration but not lysogeny in T7-like cyanophages.

Marine cyanobacteria of the genera Synechococcus and Prochlorococcus are the most abundant photosynthetic organisms on earth, spanning vast regions of the oceans and contributing significantly to global primary production. Their viruses (cyanophages) greatly influence cyanobacterial ecology and evolution. Although many cyanophage genomes have been sequenced, insight into the functional role of cyanophage genes is limited by the lack of a cyanophage genetic engineering system. Here, we describe a simple, generalizable method for genetic engineering of cyanophages from multiple families, that we named REEP for REcombination, Enrichment and PCR screening. This method enables direct investigation of key cyanophage genes, and its simplicity makes it adaptable to other ecologically relevant host-virus systems. T7-like cyanophages often carry integrase genes and attachment sites, yet exhibit lytic infection dynamics. Here, using REEP, we investigated their ability to integrate and maintain a lysogenic life cycle. We found that these cyanophages integrate into the host genome and that the integrase and attachment site are required for integration. However, stable lysogens did not form. The frequency of integration was found to be low in both lab cultures and the oceans. These findings suggest that T7-like cyanophage integration is transient and is not part of a classical lysogenic cycle.

A single-cell polony method reveals low levels of infected Prochlorococcus in oligotrophic waters despite high cyanophage abundances.

Long-term stability of picocyanobacteria in the open oceans is maintained by a balance between synchronous division and death on daily timescales. Viruses are considered a major source of microbial mortality, however, current methods to measure infection have significant methodological limitations. Here we describe a method that pairs flow-cytometric sorting with a PCR-based polony technique to simultaneously screen thousands of taxonomically resolved individual cells for intracellular virus DNA, enabling sensitive, high-throughput, and direct quantification of infection by different virus lineages. Under controlled conditions with picocyanobacteria-cyanophage models, the method detected infection throughout the lytic cycle and discriminated between varying infection levels. In North Pacific subtropical surface waters, the method revealed that only a small percentage of Prochlorococcus (0.35-1.6%) were infected, predominantly by T4-like cyanophages, and that infection oscillated 2-fold in phase with the diel cycle. This corresponds to 0.35-4.8% of Prochlorococcus mortality daily. Cyanophages were 2-4-fold more abundant than Prochlorococcus, indicating that most encounters did not result in infection and suggesting infection is mitigated via host resistance, reduced phage infectivity and inefficient adsorption. This method will enable quantification of infection for key microbial taxa across oceanic regimes and will help determine the extent that viruses shape microbial communities and ecosystem level processes.

Cyanophage host-derived genes reflect contrasting selective pressures with depth in the oxic and anoxic water column of the Eastern Tropical North Pacific.

Cyanophages encode host-derived genes that may increase their fitness. We examined the relative abundance of 18 host-derived cyanophages genes in metagenomes and viromes along depth profiles from the Eastern Tropical North Pacific Oxygen Deficient Zone (ETNP ODZ) where Prochlorococcus dominates a secondary chlorophyll maximum within the ODZ. Cyanophages at the oxic primary chlorophyll maximum encoded genes related to light and phosphate stress (psbA, psbD and pstS in T4-like and psbA in T7-like), but the proportion of cyanophage with these genes decreased with depth. The proportion of cyanophage with purine biosynthesis genes increased with depth in T4-like, but not T7-like cyanophages. No additional host-derived genes were found in deep T7-like cyanophages, suggesting that T4-like and T7-like cyanophages have different host-derived gene acquisition strategies, possibly linked to their different genome packaging mechanisms. In contrast to the ETNP, in the oxic North Atlantic T4-like cyanophages encoded psbA and pstS throughout the euphotic zone. Differences in pstS between the ETNP and the North Atlantic stations were consistent with differences in phosphate concentrations in those regimes. We suggest that the low proportion of cyanophage with psbA within the ODZ reflects the stably stratified low-light conditions occupied by their hosts, a Prochlorococcus ecotype endemic to ODZs.

Cyanobacteria and cyanophage contributions to carbon and nitrogen cycling in an oligotrophic oxygen-deficient zone.

Up to half of marine N losses occur in oxygen-deficient zones (ODZs). Organic matter flux from productive surface waters is considered a primary control on N2 production. Here we investigate the offshore Eastern Tropical North Pacific (ETNP) where a secondary chlorophyll a maximum resides within the ODZ. Rates of primary production and carbon export from the mixed layer and productivity in the primary chlorophyll a maximum were consistent with oligotrophic waters. However, sediment trap carbon and nitrogen fluxes increased between 105 and 150 m, indicating organic matter production within the ODZ. Metagenomic and metaproteomic characterization indicated that the secondary chlorophyll a maximum was attributable to the cyanobacterium Prochlorococcus, and numerous photosynthesis and carbon fixation proteins were detected. The presence of chemoautotrophic ammonia-oxidizing archaea and the nitrite oxidizer Nitrospina and detection of nitrate oxidoreductase was consistent with cyanobacterial oxygen production within the ODZ. Cyanobacteria and cyanophage were also present on large (>30 µm) particles and in sediment trap material. Particle cyanophage-to-host ratio exceeded 50, suggesting that viruses help convert cyanobacteria into sinking organic matter. Nitrate reduction and anammox proteins were detected, congruent with previously reported N2 production. We suggest that autochthonous organic matter production within the ODZ contributes to N2 production in the offshore ETNP.

Pseudo-nitzschia Challenged with Co-occurring Viral Communities Display Diverse Infection Phenotypes.

Viruses are catalysts of biogeochemical cycling, architects of microbial community structure, and terminators of phytoplankton blooms. Viral lysis of diatoms, a key group of eukaryotic phytoplankton, has the potential to impact carbon export and marine food webs. However, the impact of viruses on diatom abundance and community composition is unknown. Diatom-virus dynamics were explored by sampling every month at two coastal and estuarine locations in Washington state, USA resulting in 41 new isolates of the pennate diatom Pseudo-nitzschia and 20 environmental virus samples. We conducted a total of 820 pair-wise crosses of the Pseudo-nitzschia isolates and viral communities. Viral communities infected Pseudo-nitzschia isolates in 8% of the crosses overall and 16% of crosses when the host and viral communities were isolated from the same sample. Isolates ranged in their permissivity to infection with some isolates not infected by any viral samples and others infected by up to 10 viral communities. Isolates that were infected by the most viral communities also had the highest maximum observed viral titers (as high as 16000 infectious units ml(-1)). Titers of the viral communities were host dependent, as titers for one viral sample on eight different hosts spanned four orders of magnitude. Sequencing of the Pseudo-nitzschia Internal Transcribed Spacer 1 (ITS1) of the revealed multiple subgroups of hosts with 100% ITS1 identities that were infected by different viral communities. Indeed, we repeatedly isolated groups of isolates with identical ITS1 sequences from the same water sample that displayed different viral infection phenotypes. The interactions between Pseudo-nitzschia and the viral communities highlight the diversity of diatoms and emphasize the complexity and variability of diatom-virus dynamics in the ocean.